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We investigate the Glauber-gluon effect on the B —>■ tttt and pp decays, which is introduced via a convolution 
of a universal Glauber phase factor with transverse-momentum-dependent (TMD) meson wave functions in the 
fcr factorization theorem. For an appropriate parametrization of the Glauber phase, it is observed that a TMD 
wave function for the pion (p meson) with a weak (strong) falloff in parton transverse momentum fcr leads to 
significant (moderate) modification of the —>■ (J3° —>■ p^p°) branching ratio; the former (latter) is 

enhanced (reduced) by about a factor of 2 (15%). This observation is consistent with the dual role of the pion 
as a massless Nambu-Goldstone boson and as a qq bound state, which requires a tighter spatial distribution of 
its leading Fock state relative to higher Fock states. The agreement between the theoretical predictions and the 
data for all the B —>■ tttt and p°p° branching ratios is then improved simultaneously, and it is possible to resolve 
the B — >■ TTTT puzzle. 

PACS numbers: 13.25.Hw, 12.38.Bx, ll.lO.Hi 


I. INTRODUCTION 


The large observed -A 7r°7r° branching ratio has been known as a puzzle in two-body hadronic B meson decays, whose 
data ' [2], 

f (1.83 ± 0.21 ± 0.13) X 10"® (BABAR), 

^ 7r°7r°) = <^ (0.90 ±0.12 ±0.10) X 10-® (Belle), (1) 

I (1.17 ±0.13) X 10-® (HFAG), 


show discrepancy with the predictions obtained in the perturbative QCD (PQCD) [3] and QCD-improved factorization (QCDF) 
[4] approaches. In resolving this puzzle, one must consider the constraint from the R® —p®p® data, 

r (0.92 ± 0.32 ± 0.14) X 10"® (BABAR) , 

B{B° -A p°p°) = \ (1.02 ± 0.30 ± 0.15) x 10"® (Belle) , (2) 

I (0.97 ± 0.24) X 10-® (HFAG) , 


which, similar to the R® — 7r®7r® ones, are dominated by the color-suppressed tree amplitude C. We have carefully investigated 
the B ^ TTTT puzzle in the PQCD approach based on the kx factorization theorem [3, 5] by calculating the subleading contri¬ 
butions to the amplitude C. It was found that the next-to-leading-order (NLO) contributions from the vertex corrections, the 
quark loops and the magnetic penguin increased C, and accordingly, they increased the H® — 7r®7r® branching ratio from the 
leading-order (LO) value 0.12 x 10“® to 0.29 x 10“® [6]. At the same time, these NLO corrections increased the R® —>• p®p® 
branching ratio from the LO value 0.33 x 10“® to 0.92 x 10“® [7], which is consistent with the data in Eq. (2). Although the 
latest updates [8] of the B -A tttt analysis in the PQCD formalism have included all currently known NLO contributions, in 
particular, those to the B ^ it transition form factors [9], the agreement between the theoretical predictions and the data is still 
not satisfactory. That is, the R® -A p®p® data can be easily understood in PQCD [7] and QCDF [10], but the H® -A 7r®7r® data 
cannot. 


* Electronic address: liuxin@jsnu.edu.cn 
^ Electronic address: hnli@phys.sinica.edu.tw 
^ Electronic address: xiaozhenjun@njnu.edu.cn 

^ The latest measurement of B{B'^ —)■ = (0.90 zb 0.12 di 0.10) x 10“® with 6.7cr was released by the Belle Collaboration at ICHEP20I4 [1]. 
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The different phenomenological implication of the B ^ mr and pp data has been noticed in the viewpoint of isospin triangles 
[11], which stimulated the proposal of a new isospin amplitude with 7 = 5/2 for the latter. It has been argued [12] that the final- 
state interaction (FSI) [13, 14] could enhance the 73° —>■ 7r°7r° branching ratio through the pp —?► tttt chain. The 73° — 
branching ratio was not affected, since the tttt — pp chain is less important due to the smaller 73 tttt branching ratios. 
However, the pp —?► pp chain via the f-channel p-meson exchange was not taken into account in the above analysis. In fact, 
the p-p-p coupling is identical to the p-Tr-ir coupling in the chiral limit [15], whose inclusion will increase the 73° —p°p° 
branching ratio, and overshoot the data. Besides, the pp tttt chain is expected to enhance the 73° — tt+tt” branching ratio, 
which already saturates the data in the factorization theorems [4, 6]. Possible new physics signals from the 73 —?> tttt decays have 
been discussed in [16-18]. Similarly, a new-physics mechanism employed to resolve the 73 — tttt puzzle usually contributes to 
the B ^ pp decays, and is strongly constrained. It has been elaborated [7] that there is no satisfactory resolution in the literature: 
the subleading corrections in the factorization theorems [4, 6, 7, 10, 19] do not survive the constraints from the B ^ pp data, 
and other resolutions are data fitting, such as those by means of the jet function in the soft-collinear effective theory [20] and the 
model-dependent FSI [12-14, 21]. 

It is crucial to explore any mechanism that could lead to different color-suppressed tree amplitudes in the 73° — 7r°7r° and 
p°p° decays, and to examine whether it can resolve the 73 tttt puzzle. We have identified a new type of infrared divergence 
called the Glauber gluons [22], from higher-order corrections to the spectator diagrams in two-body hadronic 73 meson decays 
[23]. These residual divergences were observed in the kx factorization theorem for complicated inclusive processes, such as 
hadron hadroproduction [22]. They also appear in the kx factorization for 73 —Mi M 2 decays, with the M 2 meson being 
emitted from the weak vertex, which are dominated by contributions from the end-point region of meson momentum fractions. 
The all-order summation of the Glauber gluons, coupling the M 2 meson and the 73 —Mi transition form factor, generates a 
phase factor written as the expectation value of two transversely separated lightlike path-ordered Wilson lines [24]. It is noticed 
that the Glauber factor constructed in [24] is similar to the transverse-momentum broadening factor for an energetic parton 
propagating through quark-gluon plasma [25, 26]. The phase factor associated with M 2 modifies the interference between the 
spectator diagrams for C. We postulated that only the Glauber effect from a pion is significant, due to its special role as a pseudo 
Nambu-Goldstone (NG) boson and as a qq bound state simultaneously [27]. It was then demonstrated that by tuning the Glauber 
phase, the magnitude of C was increased, and the 73° — 7r°7r° branching ratio could reach 1.0 x 10“° [23]. A thorough analysis 
of 73 —Ml M 2 decays has been carried out recently, and the Glauber gluons coupling the Mi meson and the 73 —M 2 system 
were also found [28]. The resultant phase factor modifies the interference between the enhanced C and the color-allowed tree 
amplitude T. It turns out that the NLO PQCD prediction for the 73+ —?► 7r+7r° branching ratio, which receives contributions 
from both T and C, also becomes closer to the data. 

However, the Glauber phase in [23, 28] was treated as a free parameter, so it is not clear how important this phase could be. 
The postulation on the uniqueness of the pion relative to other mesons is also lacking quantitative support. According to [23], 
the Glauber phase factor is universal, depends on the transverse momenta lx of Glauber gluons, and appears in a convolution 
with decay amplitudes in the kx factorization theorem. Therefore, the universal phase factor produces different Glauber effects 
through convolutions with the distinct transverse-momentum-dependent (TMD) meson wave functions. To verify this conjecture, 
we parametrize the universal phase factor associated with the Mi and M 2 mesons as a function of the variable b conjugate to 
lx, which denotes the transverse separation between the two lightlike Wilson lines mentioned above [24]. The convolutions of 
this phase factor with the TMD pion and p meson wave functions proposed in [29], which exhibit a weaker falloff and a stronger 
falloff in the parton transverse momentum kx, respectively, indicate that the Glauber effect is indeed more significant in the 
73 TTTT decays than in the B ^ pp decays. This observation is consistent with the dual role of the pion as a massless NG 
boson and as a qq bound state, which requires a tighter spatial distribution of its leading Fock state relative to higher Fock states 
[27]. The predicted 73° 7r°7r° and 73+ — 7r+7r° branching ratios in NLO PQCD then reach 0.61 x 10“° from 0.29 x 10“° 

and 4.45 x 10“° from 3.35 x 10“°, respectively. The 73° — 7r+7r“ branching ratio decreases from 6.19 x 10“° to 5.39 x 10“°. 
Employing the same framework, we obtain the 73° —p°p° branching ratio slightly reduced from 1.06 x 10“° to 0.89 x 10“°. 
It is obvious that the agreement between the NLO PQCD predictions and the data is greatly improved for all the above modes. 

We establish the kx factorization of the 73 tttt and pp decays including the Glauber phase factors associated with the Mi 
and M 2 mesons in Sec. II. Section III contains the parametrizations of the universal Glauber phase factor, and of the intrinsic 
kx dependencies of the pion and p meson wave functions. Numerical results together with theoretical uncertainties in our 
calculations are presented. Section IV is the conclusion. 


II. FACTORIZATION FORMULAS 

In this section we derive the PQCD factorization formulas for the B{Pb) —?► Mi (Pi) M2 (P2) decay, in which the Glauber- 
gluon effect is taken into account. The B meson. Mi meson, and M 2 meson momenta are labeled by Pb, Pi, and P 2 , respec¬ 
tively, for which we choose Pb = (Pg , Pg , Og) with P^ = P^ = niBlV^, "mB being the B meson mass, and Pi (P 2 ) in the 
plus (minus) direction. The parton four-momenta k, ki, and k 2 are carried by the spectator of the B meson, by the spectator of 
the Ml meson, and by the valence quark of the M 2 meson, respectively, as labeled in Fig. 1. Specifically, we keep fc“ = xPg, 
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FIG. 1. LO spectator diagrams for the B —>■ Mi M 2 decay. 


kf = xi P^, ^2 = X 2 P 2 , and transverse components in hard kernels for 6-quark decays. For the detailed analysis of the 
Glauber divergences associated with the Mi and M 2 mesons, refer to Ref. [28]. 


A. Glauber gluons from M2 meson 

We formulate the amplitude from Fig. 1(a) for the B M 1 M 2 decay in the presence of the Glauber divergences, in which 
the hard gluon is exchanged on the right, and the Glauber gluon is exchanged on the left as shown in Fig. 2(a). The spectator 
propagator on the B meson side can be approximated by the eikonal propagator proportional to —!/(/“ + le) as I is collinear 
to P 2 , which contains an imaginary piece iTrS{l~). The propagators of the valence antiquark and quark, with the momenta 
P 2 — k 2 — k + ki — I and k 2 + I, respectively, generate poles on the opposite half-planes of as = 0. That is, the 
contour integration over l~^ does not vanish, and the Glauber gluon with the invariant mass —contributes a logarithmic 
infrared divergence f around It 0. Since a Glauber gluon is spacelike, and we are analyzing exclusive processes, 

no real gluon emissions, such as the rung gluons in the Balitsky-Fadin-Kuraev-Lipatov ladder [30], are considered. Including 
the additional Glauber divergences, we propose the Wilson links described in Fig. 3 for the modified M 2 meson wave function, 
which are motivated by the observation made in [31]: it runs from Z 2 to plus infinity along the n+ direction, along the transverse 
direction to infinity and then back to zit (the transverse coordinate of the spectator quark in the Mi meson), from plus infinity 
to minus infinity along n+ at the transverse coordinate zit, along the transverse direction to infinity and then back to the zero 
transverse coordinate, and at last back to the origin from minus infinity at the zero transverse coordinate. Moving the Wilson 
link, which runs from plus infinity to minus infinity along n+, to zit — 00 , we obtain the standard M 2 meson wave function 
[31] without the Glauber divergences. This Wilson link at the finite transverse coordinate Zit leads to the 5{l~) function. 

The modified M 2 meson wave function depends on two transverse coordinates Zi^ and Z 2 r, denoted as (ziT) Z 2 t)j where 
the dependence on 2 ^ has been suppressed. It has been shown that the Glauber gluon in the B Mi M 2 decay can be further 
factorized from the M 2 meson in the dominant kinematic region, and summed to all orders into a phase factor G{zit — zt)- We 
then have the convolution 

(ziT, Z2 t) = y d^ZTG(ziT — Zt)02(zt, Z2T + Zt), (3) 

where the definition for the two-coordinate wave function 02 (zt, Z 2 t), similar to that in [24], will be given in Eq. (8) below. 
The Wilson lines of G{zit — zt) contain the longitudinal piece, which runs from minus infinity to plus infinity in the direction 



M2 



FIG. 2. (Color online) NLO spectator diagrams for the B Mi M 2 decay that contain the Glauber divergences associated with the M 2 
meson. Other NLO diagrams with the Glauber divergences can be found in [23]. 
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FIG. 3. Wilson links for the modified wave function 4>2 ■ 


ri- = (0,1, Ot) at the transverse coordinate [24], in addition to the longitudinal piece at the transverse coordinate zit in 
Fig. 3. The above Wilson links are similar to that constructed for the jet quenching parameter in [26], which is defined as the 
average transverse momentum squared with respect to the original direction of motion that a highly energetic parton picks up, 
while traveling through a nuclear medium. If the Glauber factor contributes only a constant phase, Eq. (3) reduces to [24] 


^^(zit,Z2t) ~ exp(lS'e2) J (fzT(t>2{2T,Z2T + 2T), 
= exp(lS'e2)</'2(z2T), 


(4) 


where 4>2{^2t) denotes the standard M 2 meson wave function. The approximation in Eq. (4) with the constant Glauber phase 
has been adopted in [23, 28]. 

We route the transverse loop momentum Ij’ of the Glauber gluon through the hard gluon, the valence antiquark of the M 2 
meson, and the valence quark of the M 2 meson in Eig. 1(a). Regarding the Glauber gluon, the valence quark, and the valence 
antiquark as the partons of the M 2 meson, we assign —It, k 2 T, and —k 2 T + It to them, respectively. That is, the set of 
the Glauber gluon, the valence quark, and the valence antiquark does not carry net transverse momenta. The corresponding 
amplitude is modified into 

f (r. \2 /o /o nI [ TT^<?!>B(kT)</'l(kiT)^2(k2T,-k2T + lT)G2(lT)Tfa(kT,kiT,k2T,lT), (5) 

J {2ttY [2ttY [2ttY J {2-kY 

where the convolution in momentum fractions has been suppressed, and (fts, 4>i, and Ha denote the B meson wave function, 
the Ml meson wave function, and the hard 6-quark decay kernel, respectively. The Glauber factor G2(1 t) in momentum space 
appears as an additional convolution piece in the PQCD factorization formula for Eig. 1(a). 

The virtual gluon and the virtual quark in the hard kernel Ha have the transverse momenta Et + It — kiT and kiT — k 2 T — kT, 
respectively. We apply the variable changes kiT — It —5“ kiT and k 2 T — It —k 2 T, such that It flows through the spectator 
quark in the Mi meson, the valence quark in the Mi meson, and the valence quark in the M 2 meson. Then the It dependence 
disappears from the hard kernel, and Eq. (5) becomes 


f (Y'k.T d^kiT d^k2T 

J ( 27 r )2 ( 2^)2 ( 2^)2 


dHr 


<^B(kT)<)>l(kiT + ^t)4>2{^2T + It, 


— k2T)G2(lT)Tfa(kT, kiT, k2T)- 


(6) 


We perform the Eourier transformation of Eq. (6) by employing 


0i(kiT + lT)= J d^bi exp[z(kiT + It) ■ bi](/)i(bi), (7) 

02(k2T + It,- k2T) = J d^b2(i^b2 exp[i(k2T + It) • b^ exp[-ik2T • (-b2 - bi - b')]^2(b2,b2-F bi + b'), (8) 

G 2 ( 1 t) = J d^b'exp(ilT ■ b')exp[i5'(b')], (9) 


where the phase factor exp[i5'(b')] is a consequence of the all-order summation of Glauber gluons in 6' space [24]. Working out 
the integration over It and b 2 , and adopting 


Haihi,h2)S^^\h-hi) 


d’^'kr d^kiT d^k2T 


• b + ikiT ■ bi -f ik2T ■ ^2)Ha(kT, klT, k2T), 


( 10 ) 
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we obtain the PQCD factorization formula with the Glauber effect from the M 2 meson being included, 


d^bid^b 2 (i^b'(/)B(bi)(/<i(bi)(/) 2 (bi +b',b 2 + bi + b') exp[i5'(b')]iJa(bi, b 2 


( 11 ) 


For Fig. 1(b) with the hard gluon being exchanged on the left, we derive 


J (fhi(fh2cfh'cj)Bihi)4)i{hi)(j}2{h2 + bi + b',bi + b') exp[-i5'(b')]7Fb(bi, b 2 ). 


( 12 ) 


Note the negative phase in the factor exp[—iS'(b')], which is attributed to the Glauber gluons emitted by the valence antiquark 
of the M 2 meson. Equations (11) and (12) imply that the Glauber effect gives a strong phase to each spectator diagram for the 
color-suppressed tree amplitude. It is equivalent to route l^ through the B meson wave function in Fig. 1, under which the same 
factorization formulas will be attained. 


B. Glauber gluons from Mi meson 


We then include the Glauber gluons associated with the Mi meson, starting from Fig. 1(a). Some NLO diagrams that produce 
these types of Glauber divergences are displayed in Fig. 4. We route the transverse loop momentum 1^ of the Glauber gluon in 
Fig. 4(a) through the hard gluon, the valence antiquark of the M 2 meson, and the valence quark of the Mi meson. The spectator 
propagator on the B meson side can be approximated by the eikonal propagator proportional to —1/(/+ -f ie) as I is collinear 
to Pi, which contains an imaginary piece iTTS{l'^). The above routing of I clearly indicates that the propagators of the valence 
antiquark of M 2 and the valence quark of Mi, with the momenta P 2 — k 2 — k + ki — I and Pi — ki + I, respectively, generate 
poles on the opposite half-planes of l~ as /+ = 0. That is, the contour integration over l~ does not vanish, and the Glauber 
gluon with the invariant mass —l^ contributes a logarithmic infrared divergence as It 0. Similarly, the all-order summation 
of the Glauber divergences leads to a phase factor Gi(1t) associated with the Mi meson. The above explanation applies to the 
Glauber divergence in Fig. 4(b), and its all-order summation gives the same phase factor Gi(1t)- The reason is obvious from 
Fig. 4, where the Glauber gluon always attaches to the spectator in the B meson and the valence quark in the Mi meson. 

Assume that the Glauber gluons from the Mi meson and the M 2 meson carry the transverse momenta 1 it and I 2 T, respectively. 
The above transverse momenta are routed through the mesons, instead of through the hard kernel, so that the hard kernel has the 
same expression as in the LO PQCD approach; 1 it flows through the valence quark and the valence antiquark of the M 2 meson, 
and I 2 T flows through the Mi meson and then through the valence quark of the M 2 meson. The resultant amplitude is written as 


A,= f f ^^0^(kj,)0i(kiT + l2T -kiT - llT - I 2 T) 

“ J (27r)2 (2^)2 (2^)2 J {2Try 

X02(k2r + llT + br, —k2T — llT)Gi(liT)G2(l2T)-ffa(kr, kiT, k2T)- 


(13) 


The Fourier transformations 


Gi(1it) = J d^hei exp(iliT ■ b^i) exp [-■i5'(bsi)], 


(14) 



M2 



FIG. 4. (Color online) NLO spectator diagrams for the B Mi M 2 decay that contain the Glauber divergences associated with the Mi 
meson. Other NLO diagrams with the Glauber divergences are referred to [28]. 
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for the Glauber factor, and 


^B(kT) = J d^bBexp('ikT ■ 

4>i (kiT + hr, —kir — 1 it — br) 

= J d^b'^d^bi exp[i(kiB + hr) ■ b'J exp[-i(ki'r + br + hr) 

<p2{^2T + liT + br, —k2T — br) 

= J (fh2cfh2 exp[bk2T + br + br) ■ b^ exp[-i(k2T + br) 



(15) 

1 

1 

(16) 

(ba - h2)](t)2{K^K - b2), 

(17) 


for the meson wave functions are then inserted into Eq. (13). Note that the Glauber phases associated with Mi and M 2 differ by 
a sign for Fig. 1(a) as shown in Eqs. (9) and (14) [28]. 

We collect the exponents depending on br and br, integrate them over br and br, and obtain the 5 functions (b^i — 
b'l + bi + b 2 ) and (5(^)(bs2 + bi + b^), respectively. The next step is to perform the integration over b'^ and b^ according to 
the above 5 functions, which lead to h'^ = b^i + bi + b 2 and b^ = —bs 2 — bi. For the ky, kir, and k 2 T integrations, we still 
have Eq. (10), namely, Eb = bi. At last, we derive 


Aa= J (fhi(fh2 J(fhsi(fhs24>B{hi)^i{hsi+hi + h2,hsi+h2) 

X 02 (bs 2 +bi,bs 2 + bi +b 2 )exp [-zS'(bsi) + iS'(bs 2 )] ^^o(bi, b 2 ), 

= J (fhid'^h2 J (fhsi(fhs2^B{hi)^i{hsi+hi,hsi) 

X02(bs2,bs2 + b2)exp [-^^(bsi - b2) +iS{hs2 - bi)] Mo(bi,b2). 


(18) 


To arrive at the second expression, the variable changes bgi + b 2 —?> b^i and bs 2 + bi —> bs 2 have been employed. 

For Fig. 1(b) with the hard gluon exchanged on the left, we route br through the valence quark and the valence antiquark of 
the M 2 meson, and route br through the Mi meson and then through the valence antiquark of the M 2 meson. The resultant 
amplitude is factorized into 


Ab 


f d^k^ d^kiB d^k2T 

J (2^)2 (2^)2 (271)2 


d^liT d^bT 


(/)B(kT) 0 i(kiT + br, —kiT — br — 


br) 


x02(k2T — br, — 


k2T + br + hT)Gi{liT)G2{}2T)Hb{'kT, kir, k2T)- 


(19) 


The Fourier transformations are then applied with Eq. (17) being replaced by 

</'2(k2T — br, —k2r + br + br) 

= J (fh'2,(fh2 exp[bk2T - br) • b^ exp[-z(k2T - br - br) • (b^ - b2)]^2(b2, b^ - b2), ( 20 ) 


and we have the factorization formula 

Ab = J d'^hi(fh2 y d^b^id^b^20B(bi)^i(bsi+bi,b^i) 

X02(bs2 +b2,b^2)exp[-i5'(bsi - b2) - i5'(bs2 - bi)] iJt,(bi, ba). (21) 

The Glauber phases for the Mi meson have the same sign in Aa and At as explained before. The expressions of the hard kernels 
Ha and Hb from Figs. 1(a) and 1(b), respectively, corresponding to various tree and penguin operators, can be found in Ref. [23]. 


III. NUMERICAL ANALYSIS 


Recently, there were four works [28, 32-34] devoted to the resolution of the i? —^ tttt puzzle by enhancing the amplitude C: 

(a) In Ref. [28], Li and Mishima treated the Glauber phases as free parameters in the B ^ -ktt decays, and postulated that 
they vanish in the B ^ pp decays. When the phases associated with the Mi and M 2 mesons are both chosen as —7r/2 in 
the former, the spectator amplitudes in the NLO PQCD formalism increase, and the branching ratio becomes 

as large as 1.2 x 10“®. 
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(b) In Ref. [32], Qiao et al. significantly lowered the scale for the hard spectator interactions to the so-called optimal scale 

Qi ^ 0.75 GeV in the QCDF approach following the principle of maximum conformality, and found the 7r°7r° 

branching ratio as large as x 10“®. To justify this resolution, it is crucial to examine how the —>■ 

branching ratio is modified in the same analysis. 

(c) In Ref. [33], Chang et al. adopted large parameters pn and cj^n for the spectator amplitudes, as well as large parameters 

Pa and (f>A for the nonfactorizable annihilation ones in the QCDF framework in order to fit the tttt, -kK and KK 

data. As a consequence of the data fitting, they obtained extremely large B^ —>■ branching ratios x 10“® 

and 2.13]'[g'3g x 10“® corresponding to different scenarios. 

(d) In Ref. [34], Cheng et al. got the large color-suppressed tree amplitudes C around 0.5e“®®® and 0.6e“*®° directly 

through global fits to the data, where the former arose only from the d —?> and KK data, while the latter 

came from all the available B^^d PP data. These color-suppressed tree amplitudes resulted in the large B^ —7r°7r° 
branching ratios 1.43 ± 0.55 x 10“® and 1.88 ± 0.42 x 10“®, respectively, in the framework of flavor SU{3) symmetry. 

The experimentally observed pattern Br(i3+ —7r+7r°) > Br(i?° —tt+tt”) > Br(i7° —?► 7r°7r°) is also produced in Refs. [28, 
32]. The question on why the color-suppressed tree amplitudes are so different in the B tttt and B ^ pp decays remains to 
be answered. 

In this section we attempt to answer this question by quantitatively estimating the different Glauber effects in the B irn 
and PP decays based on the PQCD factorization formulas in Eqs. (18) and (21). As stated before, the Glauber factor is universal, 
namely, independent of the final-state hadrons, because it has been factorized from the decay processes. Nevertheless, the 
Glauber effect is not universal, as it appears through the convolution with the TMD wave functions ^i(bsi + bi,bsi) and 
02(bs2 + b 2 , bs 2 ), which possess different intrinsic b dependencies for the pion and the p meson. It will be demonstrated that 
the model wave function in [29] serves the purpose of revealing sufficiently distinct Glauber effects on the B^ and 

B'’ —>■ branching ratios. 


A. Parametrizations 


The intrinsic kx dependence of a TMD meson wave function is usually parametrized through the factor [29, 35] 

" 


k^ + m? ^ k"^ 


m 


( 22 ) 


X 1 — X 

where m = denotes the constituent quark mass, and x denotes the parton momentum fraction. Below we shall drop 

w? for simplicity. In the collinear factorization theorem one integrates a TMD wave function over kx to obtain a distribution 
amplitude. Assume that the intrinsic kx dependence appears in a Gaussian form [36], 

M'^ \ (j>M[x) 




20 


exp - 


8/3 


cjl — x) ’ 


(23) 


where Pm is a shape parameter for M = tt and p, and pM {x) denotes the standard twist-2 and twist-3 light-cone distribution 
amplitudes. Regarding the first (second) kx in Eq. (22) as the transverse momentum carried by the valence quark (antiquark) of 
the momentum fraction x (1 — x), the modified wave function is written as 

d 'k.rp d Et exp(—iky • b)^M(kTi kr), 


^3M(b',b) = 


2/32 


(27r)^ (27r)^ 

^ Pm{x) exp [-2/3^a:6'2 - 2/3|^(l - x)h'^] . 


(24) 


Our goal is to find a function 5'(b), such that the Glauber effect is large (small) for M = tt (M = p). The similar Glauber 
factor, describing the medium effect [25] in Relativistic Heavy Ion Collider physics, respects the normalization S'(O) = 0 [25]. 
If S'(b) increases with b monotonically, the real piece cos[S'(b)] takes values in both the first and second quadrants for finite b, so 
its contributions from these two quadrants cancel each other. The contribution from the third quadrant, i.e., from large b, may not 
be important due to the suppression of the exponential in Eq. (24). The imaginary piece sin[S'(b)] remains positive in the first 
and second quadrants, such that its effect always exists and becomes small only in the trivial case with S'(b) —0. Therefore, 
a monotonic function for iS'(b), which tends to enhance both the B^ tt^tt^ and branching ratios, is not preferred. A 
polynomial function or a sinusoidal function can provide an oscillatory S'(b) in b. Because the large b region is suppressed, we 
can simply parametrize >S'(b) by a sinusoidal function 


^(b) = XTT sm{pb), 


(25) 


where the tunable parameters r and p govern the magnitude and the frequency of the oscillation, and should take the same values 
for the pion and the p meson due to the universality of the Glauber factor. 
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B. Numerical results 


The following B meson wave function [3, 5] is employed in the numerical analysis, 

(I>b{x^ h) = Nbx^{ 1 — xY exp 


[ 1 

{xmB \^ 

uJ%b^ 

2 

K UJB ) 

2 


with the coefficient Nb being determined through the normalization condition 

fs 


/ dx(f)B{x, 6 = 0 ) 
Jo 


2y27Vc' 


We take the distribution amplitudes 
G/tt 


= 


2y/Wc' 

u 

2^/WZ 


[;(! — x) 


15 


1 + ^02 5(2x - 1)^ - 1 + —al 1 - 14(2a; - 1)^ + 4(2a; - 1) 


1+0 30773 - 7 ; pI ]{ 3 ( 2 a : - 1)^-1 


- o i ^ 13^3 + —p^(l + 6aJ) [> ( 3 - 30(2x - 1)^ + 35(2a^ - 1) 


2y2]v; 


(1 - 2x) 


1 + 6 ( 5773 - - 773 W 3 - —pI - -pla^ ) (1 - lOx + lOx"^) 


for the pion [37], and 


Pt^l{x) = ^^{2x-lf 

3/p 


l + ^a^'^(5(2a:- 1 ) 2-1 


l+2“2p(^5(2x-1)2-1 

3/p' 


8^6 


3fn 


^x) = ^[l + {2x-lY], rp{x) = --^(2x - 1), 


4^6 


(26) 


(27) 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


for the p meson [7, 38]. The Glauber factor is introduced only to the dominant longitudinal-polarization contribution in the 
decay. This treatment makes sense, since the Glauber effect is moderate in this mode as shown later. 

Before evaluating the i? —^ tttt and branching ratios, we explain the determination of the parameters fi-K and jjp in 
the TMD pion and p meson wave functions in Eq. (23). The parameter (3^^ around 0.40 GeV has been widely adopted in the 
literature (see for example Ref. [36]). Due to the suppression from the additional intrinsic fey dependence, we lower the shape 
parameter ujb of the B meson wave function from 0.40 GeV [ 6 ] to 0.37 GeV to maintain the NLO PQCD result for the i? —> tt 
transition form factor . The parameter /3p is not as well constrained as /3^, and we find fjp ~ P-w/i in order to maintain the 

NLO PQCD result for the B ^ p form factor Aq^^. These values of and fjp imply that the pion {p meson) wave function 
exhibits a weaker (stronger) falloff in the parton transverse momentum kx- This behavior is consistent with the dual role of 
the pion as a massless NG boson and as a qq bound state, which requires a tighter spatial distribution of its leading Fock state 
relative to higher Fock states [27]. It has been confirmed that the NLO PQCD results for all the B ^ ttit and pp decay rates are 
roughly reproduced with the above parameters, the coefficients aj = 0.115 ± 0.115, = 0.10 ± 0.10, = 0.20 ± 0.20, 

aj = —0.015, 773 = 0.015, 073 = —3, and with the chiral enhancement factor ttiq = 1.3 GeV [ 6 ], and the p 

meson decay constants fp = 0.216 GeV and /J = 0.165 GeV [39]. 

As listed in the column NLO of Table I, the NLO PQCD results for the B° 7 r° 7 r° and B'^ tt+tt^ branching ratios 
without the Glauber effect are much lower than the data, while those of the B^ —^ tt+tt” and p°p° decays overshoot the central 
values of the data. We then implement the Glauber effect, and carefully scan the r and p dependencies of the B° 7 r° 7 r° 
branching ratio. Two sets of parameters are selected, r ^ 0.47, p ~ —0.632 GeV and r ~ 0.60, p ~ 0.544 GeV, which 
give the largest B^ —>■ branching ratios 0.62 x 10“® and 0.61 x 10“®, respectively. For the former, the B^ —>■ tt+tt”, 

5+ —^ 7 r+ 7 r° and B^ —>• p°p° branching ratios are found to be 5.90 x 10“®, 3.88 x 10“®, and 1.07 x 10“®, respectively, which 
deviate from the data. For the latter, we obtain the ^ tt+tt”, —>■ and B'^ —>■ branching ratios 5.39 x 10“®, 

4.45 X 10“®, and 0.89 x 10“®, respectively, presented in the column NLOG of Table I. These outcomes show the preferred 
tendency: the B'^ tt+tt” and B'^ p°p° branching ratios decrease by 13% and 16%, respectively, while the 7 r+ 7 r° 
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TABLE I. Branching ratios from the NLO PQCD formalism in units of 10 in which NLO (NLOG) denotes the results without (with) the 
Glauber effect. 


Modes 

Data [1,2] 

NLO 

NLOG 

— >■ 7r'*’7r“ 
— >■ TT^TT® 
-K°Tr° 

B° -)■ p°p° 

5.10 ±0.19 

0 .y:O_Q 34 

0.90 ±0.16 

0.97 ±0.24 

6.19l=“^(c^B)l“:^!(aJ) 
3.351 J:?«(^b)1S:^1(oJ) 
0.29)Iq'J) (a;B)lo!o2 (®2) 
i.06t°:^®(ws)l);:l^(ag) 

5.39i):«3?(^B)l°:r5(«2) 

4.45lJ-|«(^B)l°'J®(aJ) 

0.6ll°;l«(a;B)l°;l^(aJ) 

0.89t°:?®(ws)lg:]g(a^) 


and ones increase by 33% and a factor of 2.1, respectively. The —?> tt+tt” branching ratio does not change 

much, since it is dominated by the color-allowed tree amplitude T, which is less sensitive to the Glauber effect. The ratio of 
the enhancement factor for the mode over the reduction factor for the B^ mode is about 2.5, close to 

the ratio 3 derived in Ref. [28], where the Glauber effect is assumed to be negligible in the B° —decay. Varying the 
shape parameter lub of the B meson wave function and the Gegenbauer moments ajof the pion and p meson, we estimate the 
theoretical uncertainties in our formalism given in Table I. One can see that all our predictions for the branching ratios in the 
NLO PQCD formalism with the Glauber effect match the data better. 

To have a clear idea of the Glauber effect, we present the amplitudes A{B'^ —?► and A{B° p'^) (in units of 

10“^GeV^) from Figs. 1(a) and 1(b), 


Aa.,b{B° 7r°7r°) 
Aa.b{B^ —^ p^p^) 


11.86-t9.04, -7.13 + t6.18, (NLO), 

10.80-i7.25, 7.67-i3.42, (NLOG), 

-42.44+ i24.42, 28.88-G8.07, (NLO), 

-5.78 + f4.32, -3.61-f3.23, (NLOG), 


(34) 

(35) 


respectively, associated with the four-fermion operator O 2 (they are not the full spectator amplitudes). Equation (34) indicates 
that the result of Fig. 1(a) varies a bit because of the approximate cancellation of the Glauber phases associated with the Mi 
and M 2 mesons, as shown in Eq. (18). The result of Fig. 1(b) is modified by the Glauber effect significantly with a sign flip, in 
agreement with what was found in Ref. [28]. It is obvious that the destructive interference between Figs. 1(a) and 1(b) has been 
turned into a constructive one for the B^ 7r°7r° decay. The consequence is that their sum changes from 5.53e“*°'®^ x 10“^ 
GeV^ in the NLO PQCD approach to 21.33e“*°'®^ x 10“^ GeV^ in the NLO PQCD approach with the Glauber effect. As for 
the B'^ decay, the broad distribution of the p meson wave function in b space allows cancellation to occur, which is 

attributed to the oscillation of the Glauber phase factor. This is the reason why each amplitude from Figs. 1(a) and 1(b) reduces 
as shown in Eq. (35). However, the sum of the two amplitudes does not change much relative to the change in the 5° —7r°7r° 
case. We have examined the sensitivity of the B° branching ratio to r and p, and confirmed that the predicted branching 

ratio is quite stable as long as p > 0.5 GeV, varying within only 5%. It is likely that the leading Fock state of the pion is tight 
enough to reveal the Glauber effect from the oscillatory phase factor, while other hadrons with broad spatial distributions cannot. 
We might have found plausible explanations for the dynamical origin of the Glauber phase and for the unique role of the pion 
mentioned before. 


IV. CONCLUSION 

In this paper we have performed the model estimate of the Glauber effects in the B ^ irn and pp decays in the PQCD approach 
based on the kr factorization theorem. The Glauber phase factor, arising from the factorization and all-order summation of the 
Glauber gluons for two-body hadronic B meson decays, is universal as shown in our previous work. Despite being universal, 
the Glauber factor does make distinct impacts on the B^ tAtA and B^ p^p^ branching ratios through its convolution 
with the TMD pion and p meson TMD wave functions with different intrinsic kx dependencies. It was noticed that the pion {p 
meson) wave function exhibiting a weak (strong) falloff in kx serves the purpose. These behaviors are consistent with the dual 
role of the pion as a massless NG boson and as a qq bound state, which requires a tighter spatial distribution of its leading Fock 
state relative to higher Fock states. It has been pointed out that the tight leading Fock state of the pion may be able to reveal the 
Glauber effect from the oscillatory phase factor as parametrized in Eq. (25), while other hadrons with broad spatial distributions 
cannot. 

We have demonstrated that the 7r°7r° branching ratio is enhanced by a factor of 2.1, reaching 0.61 x 10“®, while the 

BO pQpQ Qjjg j-gjiiains around 0.89 x 10“®, down by only 16%. This observation supports the fact that the Glauber effect 
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from the pion can be more significant, as postulated in [23, 28]. The tt+tt" (5+ —branching ratio is modified 

into 5.39 x 10“®, decreasing by 13% (4.45 x 10“®, increasing by 33%), such that the consistency between the NLO PQCD 
predictions and the data is improved for all the modes. The above changes are due to the facts that the Glauber phase enhances 
the color-suppressed tree amplitude by turning the destructive interference between the LO spectator diagrams into a constructive 
one, and that it also modifies the interference between the color-suppressed and color-allowed tree amplitudes. We stress that 
the B ^ TTTT puzzle must be resolved by resorting to a mechanism that can differentiate the pion from other mesons, and that 
the Glauber gluons should be one of the most crucial mechanisms. 
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